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摘 X. 为 研究 云南 伤 科 用 毒性 药材 黄 草 乌 (Radix Aconitum Vilmoriniani ) 栽培 品质 量 的 影响 因素 ， 利 用 
Illumina HiSeq 4000 高 通 量 测序 平台 对 来 自 10 个 不 同 栽培 基地 黄 草 乌 样 品 的 叶绿体 基因 组 展开 测序 , 经 过 
对 测序 数据 的 组 装 、 注 释 ， 利 用 生物 信息 学 工具 对 其 叶绿体 基因 组 特征 展开 分 析 并 构建 系统 发 育 树 。 结 果 
表明 : (1)10 个 栽培 品 的 叶绿体 基因 组 全 长 155 744—155 937 bp, 大 单 拷贝 区 和 小 单 拷贝 区 分 别 为 86 363 一 
86 548 bp、16 921—17 007 bp， 反 向 重复 区 大 小 为 26 170—26 236 bp， 均 注释 到 131 个 基因 。 (2) 序列 鉴 
定 出 60—73 个 SSR 位 点 ， 基 因 组 比较 分 析 发 现 ，10 个 栽培 品 叶绿体 基因 组 显示 出 一 定 的 扩张 ， 并 在 其 中 
发 现 了 trnK-UUU-trnQ-UUG 等 变异 热点 区 域 。 (3 ) 基 于 两 个 数据 集 的 系统 发 育 分 析 均 表明 , JS-1-4, QJ-1-2， 
LX-1-3，LJ-3-2 与 黄 草 乌 (Aconitum vilmorinianum) 亲缘 关系 较 近 ; LQ-1-3, GJ-1-3, NL-1-3, DC-2-2 和 
IRR ELO, (4. austroyunnanense) 关系 较 近 ; 基于 叶绿体 基因 组 全 序列 构建 的 系统 发 育 树 中 LJ-4-3 与 马 耳 
山 乌 头 CA. delavayi) 亲缘 关系 近 ，LJ-1-2 与 宾 川 乌 头 CA. duclouxii) 的 杀 缘 关系 较 近 ; 而 基于 蛋白质 编码 
序列 构建 的 系统 发 育 树 中 ,LJ-4-3 与 西南 乌 头 (4. episcopale) 亲 缘 关 系 近 ,LJ-1-2 M-56 ili 5535 CA. contortum? 
ne he ee de a 两 种 植 
物 ， 个 别 栽培 基地 还 掺 泥 了 乌 头 属 4conitwm) 的 其 他 物种 ， 这 可 能 是 黄 草 乌 栽 培 品 质量 不 稳定 的 因素 之 
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Abstract: To study the influencing factors on the quality of cultivated varieties of Radix Aconitum Vilmoriniani, a 
toxic medicinal material for traumatology in Yunnan Province, the chloroplast genomes of samples from ten 
different cultivation bases were sequenced by Illumina HiSeq 4000 high-throughput sequencing platform. Then 
the sequencing data were assembled and annotated followed by analysis of the characteristics of the chloroplast 
genomes by bioinformatics tools and construction of the phylogenetic tree. The results were as follows: (1) The 
full length of chloroplast genome of ten cultivated varieties from different regions was 155 744-155 937 bp, the 
large single-copy region (LSC) and small single-copy region (SSC) were 86 363-86 548 bp and 16 921-17 007 bp 
respectively, the size of the inverted repeat region (IR) was 26 170-26 236 bp. The total GC content of chloroplast 
genome of ten cultivated varieties was 38.196, with obvious AT preference, and 131 genes were annotated, 
including 85 protein-coding genes, 37 tRNA genes and 8 rRNA genes. (2) These sequences identified 60-73 SSRs 
sites, and genome comparative analysis found that the length of chloroplast genome of ten cultivated varieties 
ranged from 155 744 to 155 937 bp with moderate expansion. Some variation hot spots were found, such as 
trnK-UUU-trnQ-UUG,  trnY-GUA-trnE-UUC,  trnC-GCA-trnT-GGU,  ycfA-cemA, ycfl and ndhF, etc. (3) 
Phylogenetic analysis based on the two datasets showed that JS-1-4, QJ-1-2, LX-1-3, LJ-3-2 were closely related 
to Aconitum vilmorinianum and LQ-1-3, GJ-1-3, NL-1-3, DC-2-2 were closely related to A4. austroyunnanense. In 
the phylogenetic tree, which was constructed based on the chloroplast genomes, LJ-4-3 was closely related to A. 
delavayi and LJ-1-2 was closely related to A. duclouxii; while in the phylogenetic tree which was constructed 
based on the PCGs, LJ-4-3 was closely related to A. episcopale and LJ-1-2 was closely related to A. contortum. 
The results indicated that the hybrid origins of cultivation of Radix Aconitum Vilmoriniani included at least the 
two original plants: A. vilmorinianum and A. austroyunnanense, and other species of Aconitum even existed in 
individual cultivation bases. The mixed origins of cultivation may be one of the influencing factors on instability 

of the quality of Radix Aconitum Vilmoriniani. 
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黄 草 乌 CAconitum vilmorinianum 273-6 E fH (Ranunculaceae) 53. J& (Aconitum) & 8 k Z (Ser. Volubilia ) 
多 年 生 草本 植物 ， 分 布 于 云南 中 部 、 四 川 〈 会 理 ) 及 贵州 西部 ， 生 于 海拔 2100—2500 m 的 山地 灌 从 中 (中 
国 植物 志 编 辑 委员 会 ，1979) 。 其 药 用 部 位 为 根 ， 有 剧 毒 ， 主 要 含有 满 乌 碱 等 多 种 二 熙 类 生物 碱 ， 具 有 祛 
风 散 寒 、 活 血 止痛 、 解 毒 消 肿 等 作用 ( 李 雪 佩 等 ，2017) ， 是 云南 红 药 、 虎 力 散 和 百 宝 丹 等 云南 著名 伤 科 
类 药物 的 主要 原料 药 之 一 〈( 曾 礼 芳 等 ，2021) ， 有 着 极 高 的 药 用 价值 。 黄 草 乌 用 药 需求 大 ， 最 初 主要 依靠 
野生 资源 ， 但 过 度 采 挖 造 成 了 分 布 区 域 锐 减 和 野生 资源 枯竭 。 近 年 来 黄 草 乌 的 人 工 驯 化 栽培 获得 成 功 〈 字 
BUSES, 20120 ， 在 云南 的 沾 西 、 建 水 、 个 旧 、 东 川 、 禄 劝 、 丽 江 、 巧 家 等 地 区 有 较 大 规模 栽培 种 植 ， 有 
效 解 决 了 资源 供应 问题 。 然 而 我 们 对 几 个 栽培 基地 的 药材 进行 主要 毒性 成 分 含量 分 析 时 发 现 ， 溃 乌 碱 的 含 
量 不 稳定 ， 另 一 特征 成 分 黄 草 乌 碱 甲 的 含量 波动 范围 则 更 大 《〈 董 帅 等 ，2020) 。 黄 草 乌 作为 乌 头 类 毒性 药 
材 ， 质 量 稳定 尤为 重要 ， 因 此 其 质量 影响 因素 的 研究 很 有 必要 

黄 草 乌 药材 来 源 在 1974 年 和 1996 年 版 《云南 省 药 品 标准 》 中 规定 : 草 乌 (Radix Aconiti Vilmoriniani ) 
为 毛 草 科 植物 黄 草 乌 CA. vilmorinianum) WAS CA. austroyunnanense) 的 块根 (云南 省 卫生 局 ，1974; 
云南 省 卫生 厅 ，1996) , 2005 年 版 《云南 省 中 药材 标准 》 制 定时 仅 收 录 了 黄 草 乌 一 种 植物 来 源 ， 并 修订 药 
MHZ HRH (Radix Aconitum Vilmoriniani ) (云南 省 食品 药 中 监督 管理 局 ，2013) 。《 云 南 植 物 志 》〔 中 
国 科学 院 昆明 植物 研究 所 ,2000) 中 曾 记载 ， 黄 草 乌 与 汗 南 草 乌 在 形态 上 十 分 相近 ， 黄 草 乌 花 期 H, 
汗 南 草 乌 花 果 期 10 月 ， 块 根 均 供 药 用 ， 民 间 用 了 X ARUM. 风 竣 痪 、 跌 打 损 伤 等 疾病 ， 用 途 基 本 
相同 。 我 们 发 现 不 同 地 区 种 植 的 黄 草 乌 药材 有 黄 皮 黄 芯 、 紫 皮 黄 芯 和 紫 皮 紫 芯 几 种 情况 ， 花 的 形态 也 略 有 
不 同 ， 通 过 调研 了 解 到 ， 黄 草 乌 的 种 植 来 源 中 有 一 部 分 是 种 植 户 在 基地 附近 采 挖 的 野生 品 ， 并 未 明确 其 植 
物 基 原 。 由 此 可 见 ， 黄 草 乌 在 栽培 过 程 中 可 能 存在 种 源 混 杂 的 情况 ， 进 而 导致 药材 质量 不 稳定 。 乌 头 属 植 
物种 类 较 多 ， 种 间 界 限 模糊 ， 形 态 相似 难以 区 别 〈 朱 高 倩 等 ，2020) ， 受 外 部 环境 影响 ， 植 株 形 态 、 花 色 
花 型 、 根 的 数量 及 根 皮 颜色 等 在 种 内 也 会 有 差异 ， 仅 依靠 传统 形态 学 鉴定 、 理 化 性 质 鉴定 、 显 微 特征 鉴定 
等 方法 不 能 准确 区 别 乌 头 类 植物 CER EAS, 1986; 罗 艳 和 杨 亲 二 ,2005; 肖 培 根 等 , 2006; ARE, 2013), 
BIENNE EVI RUR EE AR E 从 源头 保障 药材 质量 

叶绿体 是 多 数 高 等 植物 和 少 部 分 藻类 及 原核 生物 进行 光合 作用 和 生命 代谢 活动 的 重要 细胞 器 ， 是 半 自 
主 性 细胞 器 ， 拥 有 自己 的 基因 组 ， 也 就 是 叶绿体 基因 组 Cchloroplast genome) 〈 王 仕 奇 等 ，2022) 。 叶 绿 
体 基 因 组 对 于 算 选 药 用 植物 密切 相关 物种 的 分 子 标 记 以 及 破译 亲缘 关系 较 近 植物 类 和 群 间 的 系统 发 育 关 系 
具有 重要 贡献 Jansen et al., 2007; Wicke et al., 2011; 姜 汶 君 等 ，2020) 。 近 年 来 新 一 代 测 序 技术 发 展 迅 速 、 
测序 成 本 逐年 降低 ， 已 经 有 越 来 越 多 的 叶绿体 基因 组 成 功 发 表 ， 推 动 着 分 子 生 物 学 研究 的 迅速 发 展 
(Shinozaki et al., 1986; Ohyama et al., 1986 ) 前 人 利用 内 转录 间隔 区 ITS( 张 富民 等 ,2003; Luo et al., 2005), 
psbA-trnH 序列 ( 朱 高 傅 等 ， 20200 等 构建 了 乌 头 属 植物 的 系统 发 育 关 系 ， 在 亚 属 水 平 上 基本 能 说 明 其 亲 
缘 关 系 ， 然 而 一 些 研究 结果 表明 亚 属 下 分 类 群 的 划分 与 传统 划分 不 一 致 ; 叶绿体 基因 组 可 以 作为 DNA 起 
级 条 形 码 对 植物 物种 进 生 了 准确 的 鉴别 ， 攀 东 昌 等 (20220 虽然 通过 叶绿体 全 基因 组 序列 构建 了 乌 头 属 的 系 


一 


© 统 进 化 树 ， 但 并 未 吉 括 马 头 属 的 大 部 分 物种 。 


本 研究 以 采集 到 的 10 个 黄 草 马 栽培 品 为 研究 材料 ， 对 他 们 的 完整 叶绿体 基因 组 序列 进行 测序 、 组 装 、 
注释 ， 揭 示 其 基本 序列 和 结构 特征 ， 同 时 结合 32 个 已 发 表 乌 头 属 物种 的 叶绿体 基因 组 序列 分 析 10 个 栽培 
品 的 系统 发 育 关系 。 拟 通过 分 析 结 果 探 讨 以 下 科学 问题 : 黄 草 乌 的 种 植 栽培 是 否 确实 存在 种 源 混 杂 的 客观 
情况 ， 进 而 导致 了 黄 草 乌 药 材质 量 的 不 稳定 。 本 下 完 将 明确 黄 章 乌 的 栽培 种 源 问题 并 为 黄 草 乌 的 基 原 植 
物 鉴定 、 规 范 化 种 植 、 合 理 开 发 利用 等 研究 提供 一 定 的 数据 。 


1 材料 与 方法 


141 实验 材料 

黄 草 乌 栽 培 品 采 自 禄 劝 县 干 海子 撤 营 盘 、 泸 西 县 三 塘 乡 、 巧 家 县 药 山 镇 、 宁 菠 县 永宁 乡 等 栽培 基地 ， 

选择 健康 植株 ， 采 和 集 新 鲜嫩 叶 并 用 干燥 硅胶 保存 作为 后 续 实验 材料 。 

K 1 10 个 栽培 品种 采集 信息 
Table 1 Collected information of ten cultivated varieties 

采集 地 编号 EN 维度 海拔 


Collection location Code ud Latitude (N) Altitude (m) 


1.2 叶绿体 DNA 提取 与 测序 
EH Bio Teke 公司 的 植物 基因 组 DNA. 提取 试剂 盒 进行 新 鲜 幼 嫩 叶 片 全 基因 组 DNA 进行 提取 ，DNA 
纯度 通过 NanoDrop 2 000 超 微 量 分 光 光 度 计 (美国 赛 默 飞 世 尔 科技 有 D 


用 于 DNA 质量 检测 , 检测 合格 后 送 往 上 海 美 吉 生物 
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HiSeq 4 000 高 通 量 测序 平台 进行 测序 。 

1.3 叶绿体 基因 组 组 装 、 注 释 和 物理 图 谱 绘制 
测序 获得 的 原始 数据 (raw data) 

数 , 除去 低 质量 的 reads 后 获得 待 组 装 序列 (clean reads) 。 参 考 乌 头 属 黄 草 乌 4. vilmorinianum (NC. 038094) 
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医药 科技 有 限 公 司 进 行 DNA. 文库 构建 ,并 采用 Illumina 


用 软件 NGS QC Toolkit (Patel & Jain, 2012) 进行 过 滤 ， 使 


JAUS 


的 叶绿体 基因 组 序列 ， 利 用 NOVOPlasty (Dierckxsens et al., 2017) 软件 对 10 个 栽培 品 的 clean reads 进行 


de novo 组 装 ，k-mer HREAN 39. XA 


在 线 注释 软件 Annotation tool GeSeq (Tillich et al., 2017) 进行 注释 ， 
并 在 Geneious R11.1.5 软件 (Kearse et al., 2012 ) 中 进行 人 工 校对 。 将 叶绿体 基因 
(Greiner et al., 2019) 在 线 工 具 ， 绘 制 叶绿体 基因 组 物理 图 谱 ， 带 有 完整 注释 信 


组 序列 导入 OGDRAW-Draw 
KT] cpDNA 序列 上 传 于 


GenBank 数据 库 中 ， 并 获得 相应 的 登录 号 (DC-2-, OP227195; GJ-1-3, OP227196; JS-1-4, OP227197; 
LJ-1-2, OP227198; LJ-3-2, OP227199; LJ-4-3, OP227200; LQ-1-3, OP227201; LX-1-3, OP227202; 
NL-1-3, OP227203; QI-1-2, OP227204) 。 


1.4 叶绿体 基因 组 特征 分 析 
更 


在 线 工 具 MISA (Beier etal., 2017) 对 各 序列 的 简单 习 
(mono-nucleotide) SSRZ10. —Z4Z Eg (di-nucleotide ) SSRZ5, —4Z HE (tri-nucleotide) SSR4, 


E 复 序列 进行 检测 ， 参 数 设 置 为 : PAER 


四 


IZER Ctetra-nucleotide) SSR Z3, AZER Cpenta-nucleotide?) SSR Z3 和 六 核 昔 酸 (hexa-nucleotide) 


SSRZ3. 
1.5 叶绿体 基因 组 比较 分 析 


分 析 。 以 黄 草 乌 为 参考 序列 ， 使 


使 用 IRscope 在 线 工 具 CAmiryousefi et al., 2018) 对 10 个 栽培 品 叶 绿 体 基 因 组 4 个 区 域 边界 进行 差异 


] mVISTA (Frazer et al., 2004) 在 线 工 具 ， 采 


] shuffle-LAGAN 模式 进行 


马 头 属 叶绿体 基因 组 全 序列 比 对 。 采 月 


H DNAsp (Rozas et al., 2017) 软件 测定 10 个 栽培 


品 叶绿体 基因 组 的 


核 苷 酸 多 样 性 Cnucleotide diversity, Pi), 参数 如 下 : 每 步 滑行 200 bpCstep-200bp) , 窗口 长 度 600 bp (window 


length-600bp) 。 
1.6 系统 发 育 分 析 


为 探讨 10 个 栽培 品 之 间 的 进化 关系 以 及 在 乌 头 属 植 物 中 的 系统 发 育 位 置 ， 从 美 


(NCBI) 下 载 了 32 个 已 发 表 的 乌 头 属 叶 绿 体 基因 组 序列 ， 此 外 ， 还 选择 了 1 种 飞 燕 草 


ERIT] 


HK 


国家 生物 信息 中 心 


E (Consolida) 植 


物 C. orientalis (NC 047292) , 1 PHE (Delphinium) 植物 D. anthriscifoliun (MK253461) 作为 外 类 


果 的 可 信和 度 ， 


和 群 ， 用 以 构建 系统 发 育 树 〈 表 2) 。 全 序列 比 对 采用 


MAFFTv.7.0 (Katoh & Standley, 2013) 进行 ， 同 时 选 
择 最 大 似 然 法 (maximum likelihood, ML) 和 贝 叶 斯 法 (bayesian inference, BI) 进行 系统 发 育 分 析 以 增加 结 


ML 树 在 软件 Phylosuite (Zhang et al., 2020) 中 进行 构建 ， 自 展 值 设置 为 1000 次 ， 贝 叶 斯 树 


同样 使 用 软件 Phylosuite 进行 构建 , 共 运 算 2 000 000 代 , 根据 马尔 科 夫 链 的 蒙特 卡 洛 模拟 算法 (markov chain 


carlo process, MCMC) ， 以 随机 树 作 为 起 始 树 ， 每 运 入 
样本 ， 且 在 平均 标准 差 降 低 至 0.01 EA FE, HAUS 
表 2 用 于 系统 发 育 树 构建 


1 000 代 
I 余 的 样本 构建 多 数 一 致 树 。 


的 已 发 表 物 种 名 录 及 其 信息 


RÉ—iX. EF burn-in 阶段 25% 的 老化 


Table 2 Published species list and its information for phylogenetic tree construction 
Genbank 序列 号 


Genus 


乌 头 属 


Aconitum 


种 中 文 名 


Chinese name of 


species 


REG 

MSk 
马 耳 山 马 头 
丽江 乌 头 
西南 乌 头 
EZ REDE 
TRU 65 3E 
德 钦 乌 头 
演 南 草 乌 
Te E 


伏 毛 铁 棒 锤 


宾 川 马 头 
显 柱 乌 头 
苍山 马 头 


PEE 
JE S 3k 
EE 
E53 


高 山 马 头 
TUE SE 


Sk 
高 帽 乌 头 


A. japonicum subsp. napiforme 


种 拉丁 名 


Latin name of species 


A. vilmorinianum 
A. hemsleyanum 
A. delavayi 
A. forrestii 
A. episcopale 
A. ramulosum 
A. bulleyanum 
A. ouvrardianum 
Aconitum austroyunnanense 
A. pendulum 
A. flavum 
A. duclouxii 
A. stylosum 
A. contortum 
A. ciliare 


A. jaluense subsp. jaluense 


A. kusnezoffii 
A. carmichaelii 
A. volubile 
A. chiisanense 
A. austrokoreense 
A. monanthum 
A. coreanum 
A. pseudolaeve 
A. finetianum 
A. longecassidatum 


A. angustius 


Genbank 
accession 
number 
NC 038094 
NC 038095 
NC 038097 
MZ959044 
NC 038096 
OM289059 
OK323949 
OM289057 
MN635745 
NC 053848 
NC 056280 
OM328070 
OM328071 
NC 038098 
NC 031420 
KT820668 
KT820670 
NC 031422 
NC 030761 
KU556690 
NC 029829 
KT820663 
NC 031423 
NC 031421 
NC 035892 
NC 036358 
NC 035894 
NC 036357 


A. reclinatum 
A. sinomontanum 
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Gymnaconitum gymnandrum 


C. orientalis 


D. anthriscifolium 


高 乌 头 
牛 扁 
ELI. i 
Consolida 
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Delphinium di 
注 : 表示 在 电子 版 《中 国 植 物 志 》 


1! 找 不 到 该 物种 的 种 


文 名 。 


MF186593 
NC 036359 
KT964698 
NC 033341 


NC 047292 


MK253461 


Note: — indicates that the Chinese species name of this species cannot be found in the electronic Flora 
Reipublicae Popularis Sinicae. 


2 结果 与 分 析 
2.1 叶绿体 全 基因 组 基本 结构 特征 和 分 类 


与 大 多 数 被 子 植 物 的 叶绿体 基因 
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(图 
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10 个 样品 叶绿体 基因 组 十 
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组 结构 相 类 似 ，10 个 栽培 品 叶绿体 基因 组 全 长 155744 〈LJ-1-2) ~ 
呈现 典型 的 四 分 体 环 状 结构 : 即 一 对 反 向 重复 区 Cinverted repeat, IR 将 
AK URS UI [X darge single copy region, LSC) F 
其 大 小 分 别 为 26 170 bp 一 26 236 bp. 85 45386 548 bp. 16 921—17 007 bp, 

呆 守 ， 不 仅 大 小 和 结构 相似 ， 而 且 GC 含量 均 为 38.1%， 都 具有 明显 的 AT 偏 
日 注释 结果 显示 〈 表 3) ，10 个 样品 叶绿体 基因 组 内 容 和 顺序 基本 一 致 ， 均 含有 
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(PCGs) 85 个 、 转 运 RNA 基因 ORNA) 37 个 和 核糖 体 RNA 基因 
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perg i 10 个 栽培 品 叶绿体 基因 组 


Chloroplast genomes of 
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二 磷酸 核 酮 糖 疲 化 酶 大 亚 基 ^ RubisCO large subunit 

光 系 统 组 装 /稳定 因素 ”Photosystem assembly/stability factors 
RNA 聚合 酶 亚 基 RNA polymerase 

核糖 体 蛋 白 小 亚 基 Ribosomal proteins (SSU) 

核糖 体 蛋 白 大 亚 基 — Ribosomal proteins (LSU) 

转运 RNAs Transfer RNAs 

核糖 体 RNAs Ribosomal RNAs 

蛋白 酶 , 成 熟 酶 clpP, matK 

其 他 基因 Other genes 

假定 叶绿体 阅读 框 Hypothetical chloroplast reading frames (ycf) 
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不 同 的 颜色 表示 不 同 功能 的 基因 。 粗 线 表 示 反 向 重复 区 (IRa 和 IRb) 的 范围 ， 它 们 将 基因 组 分 为 小 (SSC) 和 大 (LSC) 
单 拷贝 区 。 内 圈 基 因 顺 时 针 方向 转录 ， 外 圈 基 因 逆 时 针 方向 转录 。 深 灰色 内 圈 表 示 GC 含量 ， 浅 灰色 表示 AT 含量 。 

The colored bars indicate different functional group. Thick lines indicate the extent of the inverted repeat regions (IRa and IRb), 
which separate the genome into small (SSC) and large (LSC) single copy regions. Genes drawn inside the circle are transcribed 
clockwise, while those outside of the circle are transcribed counter-clockwise. The dark gray inner circle corresponds to the GC 
content, the light-gray to the AT content. 


110 个 栽培 品种 叶绿体 基因 组 物理 图 谱 
Fig. 1 Physical maps of complete chloroplast genome in ten cultivated varieties 
表 3 10 个 栽培 品种 叶绿体 基因 组 基本 信息 
Table 3 Basic information of chloroplast genome of ten cultivated varieties 
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全 基因 组 KFE N 小 单 拷贝 区 反 向 重复 区 基因 数量 
Whole LSC SSC IR Gene number 
编号 GC GC GC GC 总 基 a H rRNA tRNA 
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Total 
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( 表 4) 。 


因 , 包括 2 个 核糖 体 大 亚 基 (mp12、7p123) 、 
(yrn4.5、rrn5、rrn16、rrn23)、7 个 tRNA X 
trnR-4CG、trnV-G4C) 、1 个 NADH JB UMS (mdpB) 、5 个 未 知 功能 基因 


K 4 10 个 栽培 品种 叶绿体 


2 个 核糖 体 小 亚 基 (rps7、 rps12) 
因 GtrnA-UGC., trnI-CAU. 


基因 


组 基因 分 类 


Table 4 Functional classification of chloroplast genome genes in ten cultivated varieties 


基因 功能 
Gene 
function 
自身 复制 有 
关 
Self-replicati 
on 


基因 类 别 
Gene category 


核糖 体 蛋 白 大 亚 基 

Large subunit of ribosome 
核糖 体 蛋 白 小 亚 基 

Small subunit of ribosome 
RNA 聚合 酶 
RNA polymerase 
核糖 体 RNAs 
Ribosomal RNAs 
转运 RNAs 
Transfer RNAs 


基因 名 称 

Gene name 
rpl2(X2), rpll4, rpll6, rpl20, rpl22, rpi23(X2). 
rpl33, rpl36 
rps2, rps3, rps4, rpsT(X2), 
2), rpsl4, rpsl5. rpsl8, rps19 


rps8, rpsll, rps12(X 


rpoA, rpoB, rpoCl, rpoC2 


rrn4.5(X 2), rrnSs(X2), rrn16(X 2), rrn23(X2) 
trnA-UGC(X2), trnC-GCA, trnD-GUC, trnE-UUC, 
trnF-GAA , trnfM-CAU , trnG-GCC, trnG-UCC, 
trnH-GUG , trnl-CAU( X 2), trnI-GAU( X 2), 
trnK-UUU, trnL-CAA(X2), trnL-UAA, trnL-UAG, 


trnM-CAU, trnN-GUU( X2), trnP-UGG, trnQ-UUG, 


(96) (96) (99) (99) erof Num rRNA 
genes ber genes 
of 
Peote 
in 
codin 
g 
genes 
JS-1-4 155849 38. 86585 36.2 16924 32.6 26170 43.0 131 85 8 
QJ-1-2 155857 38.1 86582 36.2 16935 32.66 26170 43.0 131 85 8 
LX-1-3 155845 38.1 86582 36.2 16923 32.6 26170 43.0 131 85 8 
LJ-3-2 155872 38.1 86625 36.2 16921 32.5 26163 43.0 131 85 8 
LJ-4-3 155884 38. 86606 36.2 16938 032.5 26170 43.0 131 85 8 
LQ-1-3 155818 38.1 86555 36.2 17007 32.5  À 26128 43.1 131 85 8 
GJ-1-3 155818 38.1 86555 36.2 17007 32.5  À 26128 43.1 131 85 8 
NL-1-3 155811 38.1 86 548 36.2 17007 32.5 26128 43.1 131 85 8 
DC-2-2 155937 38.1 86680 36. 17003 32.5 26127 434 131 85 8 
LJ-1-2 155744 38. 86496 36.2 16944 32.66 26152 43.0 131 85 8 
10 个 栽培 品 的 基因 根据 各 功能 的 不 同 可 以 分 为 4 大 类 : OO i ERE RERH OBI EDS], @ 与 光合 作用 相关 
的 基因 ，@ 其 他 基因 ，@ 未 知 功能 基因 。 与 转录 翻译 相关 的 基因 包括 核糖 体 蛋 白 亚 基 基 因 ，RNA 聚合 酶 基 
D], rRNA 基因 和 tRNA 基因 ， 其 中 tRNA 基因 数量 最 多 ; 与 光合 作用 相关 的 基因 主要 有 光 系 统 工 基因 ， 光 
RAIA, NAD (PO H 氧化 还 原 酶 基因 ， 细 胞 色素 bf 复合 体 基 因 ，ATP 酶 基因 和 二 磷酸 核 酮 糖 羧 化 
酶 大 亚 基 基因 ， 其 中 光合 作用 系统 基因 数量 最 多 ;其 他 基因 共 6 个， 未 知 功能 基因 8 个 ， 在 这 些 基 因 中 含 
有 21 455 JL E «4 ^ rRNA 
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2.2 叶绿体 基因 组 SSR 分 析 


10 个 栽培 品种 的 叶绿体 基因 组 中 
28—37 个 单 核 背 酸 (mono-nucleotide) SSR 位 点 ，13 一 16 ^^ — EFE Cdi-nucleotide) SSR 位 点 ，6 一 
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51.5%) 。 


品 共有 的 SSR (图 2: B) 。 
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四 碱 基 重 复 单 元 AAATIATTI 和 AATG/ATTC， 五 碱 基 重 复 单元 AATAT/ATATT 是 10 个 样 
结果 表明 ，10 个 栽培 品 的 叶绿体 基因 组 中 SSRs 的 分 布 有 所 不 同 (图 2: 
区 含有 534 个 po 区 含有 88 个 SSRs, IR 区 含有 42 个 SSRs， 不 同样 品 
域 的 占 比 也 不 尽 相 同 : LSC 区 〈78.6% 一 83.1%) ，SSC [X 〈10.8% 一 14.7%) , IR 区 
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CAO 六 种 简单 重复 序列 的 数目 ，〈B) 简单 重复 序列 的 类 型 ，《〈《C) LSC, SSC 和 IR 区 域 中 SSR 的 数目 。 
(A) Number of six SSRs repeat, (B) Type of shared SSRs among the twenty-four cp genomes, (C) Number of SSRs in LSC, SSC and 


IR regions. 
图 2 10 个 栽培 品种 叶绿体 基因 组 的 简单 重复 序列 CSSR) 分 析 


Fig. 2 Analysis of simple repeat sequences (SSR) in ten cultivated varieties cp genomes 
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2.3 IR. 边界 扩张 和 收缩 分 析 

本 研究 比较 了 10 个 栽培 品 叶绿体 基因 组 的 下 区 的 边界 ， 结 果 见 图 3。10 个 样品 叶绿体 基因 组 长 度 在 
155 744 (LJ-1-2) ~155 937 bp (DC-2-2) 之 间 ，IR 区 域 的 长 度 为 26 127 (DC-2-2) 一 26 170 bp (JS-1-4， 
QJ-1-2, LX-1-3, JS-4-3) , LSC 区 域 的 长 度 为 86496 (LJ-1-2) 一 86 680 bp (DC-2-2) , SSC 区 域 的 长 
度 为 16 921 (LJ-3-2) —17 007 bp (LQ-1-3，GJ-1-3，NL-1-3) ， 所 有 样品 LSC，SSC 和 R 区 域 的 大 小 相 
fA, IR 区 与 SC 区 边界 高 度 保守 ， 但 各 样品 叶绿体 基因 组 的 IR 区 和 其 他 样品 间 仍 略 有 差异 。10 个 样品 的 
rps19 基因 均 距离 IRb [X 63 bp, KR LQ-1-3, GJ-1-3, NL-1-3, DC-2-2 4 个 样品 的 rp12 长 度 为 1 486 bp 外 ， 
其 余 样 品 的 rp12 长 度 均 为 1490 bp，ndhF 基因 距离 IRb/SSC 边界 113 一 119 bp， 所 有 样品 的 ycfl 基因 均 横 
跨 SSC 区 和 IRa X, H ycfl 基因 大 部 分 位 于 SSC 区 (4 012—4 027 bp) ， 均 有 1 286 bp 位 于 Ra 区 ; 在 
10 个 样品 中 IRA/LSC 边界 均 处 于 rpl12 和 trnH 基因 间 ，trnH 基因 距离 IRa/LSC 边界 120—126 bp. 
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3 10 个 栽培 品种 叶绿体 基因 组 LSC.. SSC, IR 边界 比较 
Fig. 3 Comparison of the junction between LSC, SSC and IR regions of chloroplast genome among the ten 
cultivated varieties 
2.4 叶绿体 基因 组 比 对 分 析 

以 黄 草 乌 为 参考 ， 使 用 mVISTA 在 线 软件 对 10 个 栽培 品 叶绿体 基因 组 进行 全 序列 比 对 ， 以 分 析 10 
个 栽培 品 叶绿体 基因 组 序列 的 差异 。 结 果 表 明 ，10 个 栽培 品 的 基因 组 高 度 保守 ， 其 中 ，IR 区 更 为 保守 ， 
编码 区 相 较 于 非 编 码 区 保守 。 此 外 ，10 个 叶绿体 基因 序列 中 存在 一 些 高 变 区 段 ， 主 要 是 基因 间隔 区 
irnK-UUU-trnQ-UUG, trnC-GCA-trnT-GGU, trnT-UGU- trnL-UAA; atpH-atpl. trnP-UGG-psaJ, rbcL-accD, 
ycfA-cemA, rpll6-rps3 以 及 基因 rp120， ycfl 等 (图 4) 。 
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上 方 的 灰色 剪 头 表 示 基 因 的 方向 ， 深 蓝 色 区 域 代表 外 显 子 ， 浅 蓝 色 区 域 代表 未 翻译 区 域 CUTRs) ， 粉 色 区 域 代表 非 编码 序 
7) CCNSO , y 轴 表 示 从 50% 到 100% 的 一 致 度 。 
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Gray arrow above the alignment indicate the direction of gene. The dark blue regions represent exons, the light-blue regions 
represent untranslated regions (UTRs), the pink regions represent Conserved Non-Coding Sequences (CNS), The y-axis represents 
the percent identity ranging from 50% to 100%. 
4 10 个 栽培 品种 叶绿体 基因 组 的 可 视 化 比 对 图 
Fig. 4 Comparison of chloroplast genomes of ten cultivated varieties 

使 用 滑动 窗口 计算 高 度 可 变 区 域 的 核 音 酸 多 样 性 ， 以 估计 所 分 析 10 个 栽培 品种 叶绿体 基因 组 中 不 同 
区 域 的 差异 水 平 。 结 果 显 示 ， 核 背 酸 多 样 性 的 变化 范围 为 0 一 0.00 915， 核 背 酸 多 样 性 大 于 0.006 的 是 基因 
间隔 区 trnY-GUA-trnE-UUC. ccsA-ndhD. petD-rpoA 以 及 基因 trnS-GCU、ycfl、rp116、ndhF， 基 因 间 隔 区 
ccsA-ndhD 变异 率 最 高 ， 核 昔 酸 多 样 性 值 为 0.00 915， 高 变异 率 位 点 主要 位 于 LSC 和 SSC 区 域 , IR 区 域 未 
发 现 变 异 率 高 的 片段 (图 5) ， 以 上 结果 与 mVISTA 软件 分 析 结 果 一 致 。 
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5 10 个 栽培 品种 叶绿体 基因 组 核 苷 酸 多 样 性 


Fig. 5 Nucleotide diversity of ten cultivated varieties chloroplast genomes 


2.5 乌 头 属 植 物 系统 发 育 分 析 

为 了 探索 不 同 种 植 基地 黄 草 乌 之 间 的 亲缘 关系 ,结合 10 个 栽培 品 叶绿体 基因 组 序列 以 及 从 NCBI 下 载 
的 32 条 乌 头 属 序列 和 2 条 外 类 群 物种 的 序列 ， 基 于 两 种 数据 集 ( 全 基因 组 序列 和 和 蛋白 质 编码 基因 序列 ) 
分 别 采用 两 种 不 同 的 方法 (ML 和 BI) 构建 系统 发 育 树 (图 6- 图 7) 。 结 果 均 显示 ，42 个 乌 头 属 物 种 聚 在 
同一 大 文 ， 整 个 乌 头 属 分 为 3 个 亚 属 ( 露 蕊 乌 头 亚 属 、 牛 遍 亚 属 和 乌 涉 亚 属 ) ， 露 芍 乌 头 亚 属 为 基部 类 和 群 ， 
10 个 栽培 品 均 属 乌 头 亚 属 。 以 不 同 数据 集 建 立 的 进化 树 在 拓 朴 结构 上 存在 细微 差异 ， 基 于 叶绿体 全 基因 组 
构建 的 系统 发 育 树 均 显示 ，LJ-1-2 与 铁 棒 锤 CA. pendulum) 等 3 个 物种 聚 为 分 支 as， 支 持 率 (99/1) ,与 
宾 川 乌 头 (4. duclouxii) 亲 缘 关 系 较 近 ; LQ-1-3、GJ-1-3、NL-1-3 和 DC-2-2 与 演 南 草 乌 (4. austroyunnanense ) 
聚 为 分 支 bp， 支 持 率 (100/1) ; LJ-4-3 与 马 耳 山 乌 涉 (A. delavayi) ~ WSA CA. forrestii) 等 3 个 物种 
以 (100/1) 的 支持 率 聚 为 分 支 c， 与 马 耳 山 乌 头 CA. delavayi) 亲缘 关系 较 近 ; JS-1-4、QJ-1-2、LX-1-3 以 
及 LJ-3-2 与 黄 草 乌 CA. vilmorinianum) 和 瓜 时 乌 头 〈4. hemsleyanum) 聚 为 分 支 4， 与 黄 草 乌有 着 较 近 的 
亲缘 关系 ;而 基于 和 蛋白质 编 码 基 因 构 建 的 系统 发 育 树 均 显示 ，LJ-1-2 5Eitiili 3k CA. contortum). 聚 为 分 支 
a, XR (88/1) ; LQ-1-3、GJ-1-3、NL-1-3 和 DC-2-2 与 滇 南 草 乌 CA. austroyunnanense) 聚 为 分 支 b， 
支持 率 (10001) ; LJ-4-3 与 西南 乌 头 CA.episcopale) 、 马 耳 山 乌 头 CA. delavayi) 等 4 个 物种 以 (74/1) 
的 支持 率 聚 为 分 支 c， 与 西南 乌 头 亲缘 关系 较 近 ;JS-1-4、QJ1-2、LX-1-3 以 及 LJ-3-2 与 黄 草 乌 〈4. 
vilmorinianum) 聚 为 分 文 d。 
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If the species name is not marked in Chinese, it means that the Chinese name of the species cannot be queried in the electronic Flora 


Reipublicae Popularis Sinicae. The same below. 


6 基于 44 个 物种 叶绿体 全 基因 组 序列 构建 的 最 大 似 然 MLO 和 贝 叶 斯 CBD 系统 进化 树 
Fig. 6 Phylogenetic tree of 44 species using maximum likelihood (ML) and bayesian inference (BI) methods 


based on chloroplast genome sequences 
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7 基于 44 个 物种 蛋白 质 编码 基因 序列 构建 的 最 大 似 然 (ML) 和 贝 叶 斯 BDO 系统 进化 树 
Fig. 7 Phylogenetic tree of 44 species using maximum likelihood (ML) and bayesian inference (BI) methods 


based on protein-coding gene sequences 
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3 讨论 


药材 质量 的 影响 因素 有 内 因 和 外 因 ， 主 要 是 植物 基 原 和 环境 因素 两 方面 。 黄 草 乌 作为 云 药 伤 科 要 药 ， 
同时 又 是 剧 毒 类 药材 ， 种 植 基 原 的 混乱 可 能 会 造成 化 学 成 分 种 类 和 含量 的 差异 ， 进 而 影响 药材 质量 的 稳定 
性 ， 用 药 安全 隐患 极 大 。 目 前 有 学 者 为 了 黄 草 乌 的 合理 安全 利用 ， 对 其 化 学 成 分 〈 汪 焕 芹 等 ，2014) 、 
理 作用 ( 温 玉 莹 等 ，2019)〉 和 生 药 学 鉴定 ( 王 丽 等 ，2020) 5 等 方面 展开 了 研究 ， RILA RE E G JUR m 
基 原 植物 叶绿体 基因 组 的 报道 。 本 研究 利用 高 通 量 测序 技术 和 生物 信息 学 分 析 方 法 对 10 个 栽培 品 的 叶 绿 
体 全 基因 组 特征 和 系统 发 育 关系 展开 了 研究 探讨 ， 结 果 显 示 ，10 个 栽培 品 的 叶绿体 基因 组 大 小 相似 ， 均 以 
典型 的 、 高 度 保守 的 四 分 体 结构 存在 ， 包 括 四 个 部 分 ，1 个 LSC 区 ，1 个 SSC 区 和 1 对 IR 区， 叶绿体 基 
因 组 长 度 差异 仅 为 190 bp, GC 含量 均 为 38.1%， 具 有 明显 AT 偏好 性 ， 均 注释 到 131 个 基因 ， 这 与 乌 头 属 
其 他 物种 的 叶绿体 基因 组 研究 结果 十 分 相似 CLiu et al., 2020; Zhang et al., 2021) ， 说 明 叶绿体 基因 组 进化 
速率 较 慢 ， 在 同属 植物 间 表 现 出 相对 保守 性 。 

人 多 态 性 高 、 易 复制 且 具 有 母系 遗传 模式 
等 优点 而 被 广泛 应 用 于 遗传 多 样 性 分 析 和 分 子 标记 辅助 育种 等 方面 的 研究 (Pauwels et al., 2012) 。 本 研究 
通过 在 线 分 析 软 件 在 10 个 栽培 品 的 叶绿体 基因 组 中 夫 URSLGU-73 ^- SR As, KOERREN 
最 大 ， 其 次 为 二 核 苷 酸 和 三 en 
AT/TA 是 二 核 苷 酸 主 要 重复 类 ， 三 核 背 酸 主 要 重复 类 型 则 是 AATATT， 这 与 李 雪 佩 等 (2018) 研究 的 云 
南 省 武 定 县 的 黄 草 icol: d nd c: 
提供 候选 分 子 标 记 。 虽 然 叶 绿 体 基 因 组 高 度 保守 , 但 是 基因 结构 和 大 小 在 10 个 栽培 品 叶绿体 基因 组 的 IR/SC 
边界 也 有 一 定 不 同 ， 对 10 个 栽培 品 进行 全 序列 对 比 ， 发 现 他 们 的 LSC，SSC IIR 区 域 的 大 小 存在 细微 差 
t, 分别 为 86 496—86 680 bp，16 921—17 007 bp, 26 127—26 170 bp， 显 示 出 适度 的 扩张 ， 相 对 于 LSC 


im 


区 和 SSC X, IR 区 更 为 保守 ， 这 可 能 是 由 于 IR. 区 的 重复 基因 会 在 突变 发 生 时 ， 通 过 基因 转换 对 其 进行 修 
正 ， 因 此 减弱 了 此 区 变异 发 生 的 频率 (Chen etal., 2015) ， 所 有 样品 rpsl9. ndhF. ycfl 和 trnH 基因 的 大 


小 和 分 布 在 四 个 边界 间 的 差异 较 小 。 基 于 mVISTA 和 DNAsp 软件 分 析 结 果 表 明 , 在 10 个 样品 存在 一 些 变 
异 率 较 高 的 区 段 , 如 trnK-UUU-trnQ-UUG., trnY-GUA-trnE-UUC, trnC-GCA-trnT-GGU, trnT-UGU-trnL-UAA, 
atpH-atpl, trnP-UGG-psaJ, rbcL-accD, ycfA-cemA, rpll6-rps3 以 及 基因 rp120, ycfl, ndhF 等 , Park 等 (2017) 
和 Meng 等 (2018) 也 有 相似 的 研究 结果 。 在 以 往 的 研究 中 ,ndhF-trnL, psbA-trnH, psbD-trnT 和 trnT-trnL 
等 片段 被 用 于 乌 头 属 的 系统 发 育 研究 (Hong etal., 2017) ， 但 乌 头 属 种 间 亲 缘 关 系 未 能 得 到 很 好 的 解决 ， 
今后 这 些 变异 热点 区 段 ， 或 许可 以 作为 新 的 DNA. 条 形 码 ， 用 于 黄 草 乌 及 乌 头 属 植物 的 物种 鉴定 和 系统 发 
育 关系 等 方面 的 研究 。 
与 已 有 的 条 形 码 片段 相 比 ， 叶 绿 体 全 基因 组 序列 具有 更 强 的 分 辩 率 ， 可 以 阐明 较 低 阶 元 的 系统 发 育 关 

系 (Nie etal., 2012) 。 本 研究 基于 两 种 数据 集 采 用 不 同 的 方法 构建 的 系统 发 育 树 均 表明 ， 乌 头 属 和 2 个 种 
分 为 三 个 大 的 分 支 〈 乌 头 亚 属 、 牛 扁 亚 属 和 露营 乌 头 亚 属 ) ， 乌 头 亚 属 与 牛 扁 亚 属 构成 姊妹 类 群 ， 露 划 乌 
FI 头 亚 属 为 乌 头 属 的 基部 类 群 ，10 个 栽培 品 均 属 乌 头 亚 属 ， 但 仅 有 JS-1-4.. QI-1-2.. LX-1-3 和 LJ-3-2 4 4-3 
= 培 种 与 黄 草 乌 亲 缘 关 系 较 近 ，LQ-1-3、GJ-1-3、NL-1-3 和 DC-2-2 与 滇 南 草 乌 亲 缘 关 系 近 ， 丽 江 的 LJ-4-3 
在 叶绿体 全 基因 组 序列 构建 的 发 育 树 中 与 兴安 乌 头 系 (Ser. Ambigua) 的 马 耳 山 乌 头 亲缘 关系 更 近 ， 在 蛋白 

质 编 码 基因 序列 构建 的 发 育 树 中 则 与 草 乌 头 系 〈Ser. Volubilia) 的 西南 乌 头 关系 近 ; 丽江 另 1 个 样品 LJ-1-2 

在 叶绿体 全 基因 组 序列 构建 的 发 育 树 中 与 保山 乌 头 系 〈Ser Bullatifolia) 中 的 宾 川 乌 头 关系 更 近 , 在 蛋白 质 

编码 基因 序列 构建 的 发 育 树 中 则 与 显 柱 乌 头 系 〈Ser. Stylosa). 的 苍山 乌 头 关 系 近 ， 这 说 明 非 编码 区 在 植物 

系统 进化 中 可 能 也 起 到 了 一 定 的 作用 。 可 见 ， 除 丽江 地 区 的 栽培 品 来 源 较为 混乱 外 ， 其 他 地 区 的 栽培 种 从 

亲缘 关系 上 基本 分 为 黄 草 乌 和 汗 南 草 乌 两 大 类 ， 这 说 明 目 前 各 基地 种 植 的 黄 草 乌 主要 存在 两 种 植物 来 源 。 

此 外 ， 研 究 结果 还 显示 蔓 乌 头 系 、 兴 安 乌 头 系 、 显 柱 乌 头 系 、 保 山 乌 头 系 和 短 柄 乌 头 系 等 的 部 分 类 群 交 错 
在 一 起 ， 这 与 中 国 植物 志 中 乌 头 属 系 的 划分 并 不 一 致 ， 前 人 《“〈 张 富民 等 ，2003) 也 得 出 了 相似 的 结果 ， 这 
提示 依靠 植物 形态 的 经 典 分 类 系统 可 能 存在 一 些 问题 。 鉴 于 乌 头 属 植物 在 生长 过 程 中 性 状 的 多 变性 ， 仅 依 
靠 表 型 特征 对 其 进行 区 别 分 类 显然 已 不 足够 ， 这 也 可 能 是 黄 草 乌 在 引种 驯化 的 过 程 中 出 现 种 源 混杂 情况 的 
原因 。 虽 然 目 前 乌 头 属 大 部 分 植物 已 经 完成 叶绿体 全 基因 组 测序 ， 但 有 关 该 属 植物 的 系统 发 育 学 研究 还 比 
较 少 ， 本 研究 推断 的 系统 关系 能 够 为 乌 头 属 植 物 的 物种 鉴定 和 系统 进化 等 方面 的 研究 提供 重要 信息 ， 但 乌 
头 属 植物 种 类 繁多 ， 仍 然 需要 更 多 该 属 植物 加 入 到 系统 发 育 的 研究 中 且 结 合 更 多 技术 手段 ， 才 能 得 到 更 全 
、 更 合乎 实际 的 进化 谱 。 
综 上 所 述 ， 本 研究 从 叶绿体 基因 组 的 角度 认为 ， 目 前 云南 省 各 大 栽培 基地 所 种 植 的 黄 草 乌 确实 存在 种 

源 混杂 的 客观 现象 ， 主 要 有 黄 草 乌 和 溃 南 草 乌 两 种 植物 ， 但 现行 质量 标准 仅 有 黄 草 乌 一 种 植物 来 源 ， 因 此 
对 于 黄 草 乌 栽培 药材 的 质量 问题 , 首先 应 规范 其 种 源 , 而 对 于 溃 南 草 乌 能 否 也 应 作为 黄 草 乌 药材 来 源 收 载 ， 


通过 研究 他 们 化 学 成 分 和 药理 活性 等 方面 的 异同 综合 评判 。 
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